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ABSTRACT 

The origin of tilted disks in cataclysmic variables is explained in terms of a model involving the 
stream-disk interactions. Tilted, precessing disk causes periodically variable asymmetry in the irra- 
diation of the two hemispheres of the secondary component, resulting in variable vertical component 
of the velocity of the stream. The following stream-disk interactions provide additional vertical ac- 
celeration to disk elements needed to produce and maintain disk tilt. Predictions based on this model 
compare favorably with observations. 
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1. Introduction 

Negative superhumps are quite common among cataclysmic variables. They 
are present in dwarf novae of the SU UMa type during their superoutbursts and also 
among systems with stationary accretion - the permanent negative superhumpers 
(see Patterson 1999, Wood and Burke 2007, Olech et al. 2009 and references 
therein). Their periods are shorter than the orbital periods, the corresponding "pe- 
riod deficit" £ being correlated with the orbital period. There are examples of 
negative superhumpers showing also the common superhumps, the two types either 
being present simultaneously or switching from one type to another. In those cases 
the "period excess" of common superhumps and the "period deficit" of negative 
superhumps are correlated (cf. Olech et al. 2009). 

The commonly accepted interpretation of negative superhumps explains them 
in terms of a tilted, precessing disk (cf. Patterson et al. 1993, Harvey et al. 1995, 
Patterson 1999, and references therein). Supporting this interpretation are: (1) the 
presence of light variations with Pprec observed commonly in such systems and (2) 
the fact that precession periods obtained from the observed values of and Porb 
agree reasonably well with theoretical predictions (see Larwood et al. 1996 and 
references therein). 
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Patterson et al. (1997) were the first to suggest that negative superhumps could 
be due to the "spot" produced by the variable stream impact as it transits across 
the surface of a tilted disk. More recently this was confirmed by Wood and Burke 
(2007) who used the 3D SPH simulations of a tilted, precessing disk to produce 
light curves closely resembling the observed superhumps (see also Montgomery 
2009 and Wood et al. 2009). 

The origin of the disk tilt, however, continues to be unknown. In particular, 
models proposed for other, more exotic objects (such as HZ Her or SS 433), are not 
applicable to the case of disks in cataclysmic variables (cf. Wood and Burke 2007 
and references therein). 

In the present paper we propose to explain the origin of tilted disks in cat- 
aclysmic variables in terms of a model involving the stream-disk interactions in 
a situation when the tilted, precessing disk causes variable irradiation of the two 
hemispheres of the secondary component and the resulting stream has periodically 
variable vertical velocity component. 

We begin, in Section 2, with definitions and formulae to be used in further 
sections. Section 3 describes the irradiation controlled mass outflow and the next 
two sections are devoted to a detailed discussion of the resulting stream: its tra- 
jectory (Section 4) and its collision with the surface and edge of the disk (Section 
5). The model, involving the stream-disk interactions, is presented in Section 6, its 
predictions being compared with observations in Section 7. 



2. Definitions and Formulae 



The precession period of a tilted disk is related to the orbital period and the 
negative superhump period by 

' ' ' (!) 



Pprec fnsh ^orb 

Accordingly we have 



^prec — ^nsh ^orh i (2) 

with the zero-points of phases being defined as follows: ^orb = at conjuction (in 
particular - at eclipse) and ^nsh = at superhump maximum. 
It will also be useful to recall that 



where 



d^orb 1 - e ' 



8 = ^"'•^ ^"^^ , (4) 
Porb 
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is tiie negative superhump period deficit. 

Let us now consider a disk with radius r^, geometrical thickness z/r, and tilt 
angle 6. Its geometry (see Fig.l) is described by the following set of equations: 

Zd = - Zo cos9 , (5) 
where 6 is the position angle on the surface of the disk and 

Zo = rd sin 5 , (6) 

and 

= Zrf + Azrf , Zb = Zd - ^ , (7) 




Fig. 1. Geometry of the tilted disk with 5 = 3° (see Section 7.1) and z/r = 0.10 at 6 = 0. Dotted 
Une marks the orbital plane. 

Two comments are worth adding here. First, that the specific value of the disk 
radius is actually irrelevant since in what follows we will be dealing only with 
angles. Secondly, that all angles will be expressed here in units of phase, so that 
6 = 1 corresponds to 2% or 360° . 

Let us now consider an element of the outer disk. Its motion in the z-coordinate 
is described by Eq.(5) with 

2% 

e = ^<, (9) 

where Pj is the period of oscillations (equal also to the period of revolution). 
The acceleration, calculated from Eqs.(5) and (9), is 



(fz f2nV „ ^ dzo f2n\ . „ 

The second term on the right hand side of this equation implies that to pro- 
duce and maintain disk tilt we need a mechanism capable of providing additional 
acceleration of the form 



~ sin6 . 



(11) 
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3. 



The Irradiation Controlled Mass Outflow 



It has recently been shown (Smak 2008) that superoutbursts of Z Cha (and 
most hkely of other dwarf novae) are due to a major enhancement in the mass 
outflow rate. This provides a strong argument in favor of the concept of irradia- 
tion controUed mass outflow. The details of this mechanism, however, still remain 
controversial (see Smak 2009b). The problem here is connected with the fact that 
the equatorial parts of the secondary, including L i , are in the shadow cast by the 
disk. The crucial question then is whether the material from irradiated regions 
flowing towards L i is still hot enough after reaching this point to produce substan- 
tial modulation of the mass outflow rate. Assuming that this is indeed the case it 
was possible to propose a new interpretation for superhumps (Smak 2009b). The 
simplicity and self-consistency of this interpretation could, in fact, be treated as an 
indirect argument in favor of this assumption. 

In what follows we also assume that the outflow from L i is controlled by ir- 
radiation and, in particular, that in the case of variable irradiation there is a time 
delay between irradiation and the resulting dissipation at the point of impact (see 
Smak 2009b): 



where ^fiow is the time needed for the flow from irradiated regions to reach Li, 
and i!Ustr is the time needed for the stream to reach the point of impact. The flow 
time ^^fiow depends, in general, on several parameters; in particular it depends on 
the distance between L i and the shadow boundary. 



When the disk is coplanar with the orbital plane the initial velocity vector of 
the stream at L i is of the form Vo = [v;c,o , Vy^o , 0] . Let us, however, consider the 
case of a tilted disk when its part facing the secondary is tilted below the orbital 
plane (in the notation used in Section 6 this corresponds to %2,d = 0)- Compared 
to the coplanar case, the shadow boundary on the top hemisphere of the secondary 
is now closer to Li , while that on the bottom hemisphere - further away from Li . 
As a result, the contribution to the mass outflow from the material flowing from the 
top hemisphere is dominant and, consequently, the initial velocity vector now has 
anon-zero vertical component: v^^o < 0. 

To study the consequences of such a situation we calculate stream trajecto- 
ries in three dimensions. For the initial stream velocity at Li we adopt vari- 
ous combinations of its components (in dimensionless units): Vx^o = 0.01 —0.05, 
Vy^o = (-0.01) - (-0.05), and v^,o = (-0.02) - (-0.05). Shown in Fig.2 are re- 
sults obtained for the mass ratio q = 0.3 which corresponds to the typical periods 



A? = Atfiow + A?^ 



str 1 



(12) 



4. Stream Trajectories 
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of negative superhumpers, as for example those listed below in Table 1. (Results 
obtained with other mass ratios were qualitatively the same). 
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Fig. 2. (a) Stream trajectories in projection on the x-y plane, (h) Stream trajectories in projection on 
the x-z plane, (c) Vertical velocity component of the stream, (d) The "flight" time Atgir in units of 
the orbital period. Vertical dotted line in all plots marks the location of the disk edge. 

There are three important conclusions to be drawn from Fig. 2: (1) The stream 
remains very close to the orbital plane. (2) The "flight" time from Li to r = is 
Atstr ~ 0. ISPorb and depends only weakly on the initial conditions. (3) At the point 
of impact the sign of the vertical velocity component of the stream is opposite to 
that of the initial velocity. 

5. The Stream Overflow 

The stream overflow was originally expected (cf. Smak 1985, Hessman 1999) 
to occur when the disk is geometrically thin, i.e. mainly in quiescent dwarf novae. 
Evidence is now available, however, for a substantial stream overflow in Z Cha 
during its superoutbursts (Smak 2007, 2009a). Besides, there is also theoretical 
evidence (Kunze et al. 2001) suggesting that this is indeed a much more common 
phenomenon. 
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Using disk geometry described in Section 2 we can calculate the fraction ft of 
the stream material overflowing the top part of the disk and the fraction fg colliding 
with disk edge. We adopt: z/r = 0.10 (cf. Smak 1992), 6 = 3° (see Section 7.1) 
and the density distribution in the stream given by the two-dimensional Gaussian 
formula with a/Azj = 0.5, 1.0, and 2.0. With these assumptions we have 

ft,e = / ^ exp (-zV2a2) dz , (13) 

Zl 

where the integration hmits are: zi = Zt and Z2 = °° for ft , and z\ = Zb and zi = Zt 
for fe. 
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Fig. 3. The fraction // of the stream material overflowing the top part of the disk (broken Unes) and 
the fraction colliding with its edge (solid lines) are shown as a function of the position angle for 
5 = 3° and for three values of o . 

Results are presented in Fig.3. As could be expected, substantial overflow oc- 
curs when a is comparable to Az^ . In particular, f becomes larger than fe for 
a > 2diZd ■ Furthermore, variations of /( show only one maximum (and one mini- 
mum) per cycle, while those of /g - two maxima (and two minima) per cycle (we 
shall return to this in Section 7.2). 

Additional calculations were also made to reproduce the simple situation con- 
sidered by Wood and Burke (2007) with the stream fully overflowing one side of 
the disk (f ^ I). We find that this would require either much larger values of a 
(which is rather unlikely) or tilt angles much larger than 5 ~ 6° (in fact. Wood and 
Burke used 6 = 5°, but their disk was much thiimer). We shall return to such a 
situation in Section 7.1. 

6. The Model 

The stream, colliding with the disk with non-zero vertical velocity component, 
is an obvious source of extra acceleration of disk elements. We propose that, under 
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suitable conditions, this can be the mechanism capable of producing and maintain- 
ing disk tilt. 

The geometry of irradiation depends on the relative orientation of the tilted disk 
with respect to the secondary component. The position angle of the secondary with 
respect to the observer is 

Q2,obs = ¥orb ■ (14) 

The asterisk * is used here (and below) to emphasize that the phase refers to the 
moment of irradiation of the secondary component. The position angle of the low- 
est point of the disk, corresponding to 9 = (as defined above), with respect to the 
observer is 



Qd,obs = <t>prec " , (15) 

where 9o is the precession phase at which the point with 9 = is facing the ob- 
server. Note that Bd,ohs , which changes due to retrograde precession, is counted in 
the direction opposite to Qi^obs ■ With this definition, the position angle of the disk 
point with 9 = with respect to secondary component is 

Qd,2 = Qd,obs + Q2,obs = {<^*prec + Krb) " = <^lsh " Qo ■ (16) 

Consequently, the vertical component of the terminal stream velocity at the 
point of impact v,mp can be expressed as a function of (|)*^^ . Assuming, for sim- 
plicity, that its variations are cosinusoidal we can write 

Vimp = Vimp,o COS ( (|)*^;, - 9o ) , (17) 

where vtmp.o > and corresponds to (|)*^,, = 9o (see Fig.2c and Eq.l6). The addi- 
tional acceleration produced by the stream is then 



az,str = M Vimp,o COS ( ^^.^^ - 9o ) ~ sin [ ( C,;, - 9o ) + 0.25 ] . (18) 

The moment of impact is delayed with respect to the moment of irradiation by 

Af (see Eq.l2), the corresponding phase delays being ^(Sforh = /Porb (see Fig.2d) 
and A(|)„^/, = A(|)or&/(l — e) (see Eq.3). The phase of impact is then 

<^nsh = ¥nsh + ^<9nsh = + ^ , (19) 

and the position angle on the disk at the point of impact 

6 = Ksh - 60 + A^nsh + a = rnsh -Qo + ^+ a, (20) 

where a is the position angle of the impact point with respect to the line joining the 
two components. For q = 0.3 and at r = we have a 0.014 (in phase units). 
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The acceleration needed to produce and maintain disk tilt (Eq. 11 in Section 2) 
must then be 



a, ~ sm 



( - 9o ) + y— ^ + a 



(21) 



Comparing Eqs.(18) and (21) we immediately conclude that our proposed mech- 
anism is most efficient when 



1-8 



+ a = 0.25 , 



(22) 



or, using 8 0.02 and a 0.014, when 



K<^orb ~ 0.23 . (23) 

In general the efficiency of this mechanism is described by the sign and value 
of the integral 



1 

W = Jcos (C;, - Oo) sin 



nsh 



1 



71 sm I h a 



V 1-8 

which implies that the mechanism is effective (W > 0) for 



(24) 



- 0.02 < A^orh < 0.48 . (25) 

The value of A(^orh , estimated earlier (Smak 2009b) for three systems (lY UMa, 
DV UMa and OY Car) showing common superhumps, was found to be A^orb = 
^/Porb ~ 0.66 which is clearly outside the range required by Eq.(25). This is not 
surprising, however, since those three objects do not show negative superhumps. 

In the case of negative superhumpers with tilted disks the shadow boundary 
on the top (or bottom) hemisphere of the secondary is much closer to Li and 
this makes the flow time (Eq.l2) much shorter. Consequently the value of Ai^orb 
can easily be small enough to fulfill condition imposed by Eq.(25). To illustrate 
this point let us consider the extreme case, when the shadow boundary is close to 
Li (see Section 7.1). In such a case At flow and A^orb becomes as small as 
Atstr/Porb ~ 0-15 (scc Eq.l2 and Fig.2d). 



7. Model Predictions 



7.1. Disk Tilts 

Let us consider the situation when the tilt is so large that the vicinity of L i 
is fully exposed to direct irradiation. In such a case the vertical component of 
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the stream velocity becomes Vj-^o = and the stream mechanism discussed in the 
previous Section no longer operates. This means that there is a natural upper Umit 
to the tilt angle given by dmax = arcsin(z/r) . Using values of z/r = 0.10 — 0.12 
typical for hot disks we predict 5,„cx ~ 6 — 7° . 

To test this prediction we analyze the best documented examples of negative 
superhumpers which show periodic light variations with Pprec ■ They are listed in 
Table 1, where the second and third columns contain the semi-amplitudes A of 
those variations and the orbital inclinations; in the case of non-eclipsing systems, 
with no estimates of inclinations being available, we adopt a range: / = 20 — 60° . 



Table 1 

Semi- Amplitudes and Tilt Angles 



Star 


A (mag) 


j(deg) 


Refs. 


6(deg) 


PXAnd 


0.20 


74 


1,2 


3 


V603 Aql 


0.07 


20 


3 


7 


TT Ari 


0.07 


29 


4,5,6,7 


4 


TV Col 


0.20 


70 


8,9,10 


3 


V751 Cyg 


0.05 


20-60 


11 


5-1 


V1084 Her 


0.08 


20-60 


12 


6-2 


V442 Oph 


0.10 


67 


12,13 


2 



References: 1. Stanishev et al. (2002, Fig.5), 2. Thorstensen et al. (1991), 3. 
Patterson et al. (1997, Table 1), 4. Semeniuk et al. (1987, Fig.9), 5. Kraicheva 
et al. (1997), 6. Kraicheva et al. (1999), 7. Wu et al. (2002), 8. Barrett et al. 
(1988, Fig.4), 9. HelUer (1993), 10. Retter et al. (2003, Table 4), 11. Patterson et 
al. (2001, Fig.3), 12. Patterson et al. (2002, Figs.8 and 10), 13. Ritter and Kolb 
(1998). 

The semi-amplitude A (in magnitudes) can be written as 

A = ^8, (26) 

di 

where dMy/di describes the dependence of the observed luminosity of the disk on 
its inclination (or - generally - on the viewing angle). 

At inclinations i < 75° it is sufficient to use the flat disk approximation giving 

Ld{i) = <Ld> { \ — u + u COS? ) COS? , (27) 
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where u is the limb darkening coefficient for which we adopt m = 0.6 . Turning to 
magnitudes we get 

dMv , (1— m + 2m cos / ) sinj\ . ^ 

— ^ = 0.01895 \- ^ mag/deg . (28) 

ai ( 1 — M + M cos; ) COS? 

The resulting values of 8 are listed in the last column of Table 1 . As one can 
see they are smaller that b^ax predicted above, the typical tilt being 5 3° . 

7.2. The Superhump Light Curves 

The negative superhumps are thought to be due - primarily - to the collision 
of the stream with the surface of the tilted, precessing disk. There are two other 
effects, however, which must be considered. First, that part of the stream collides 
with disk edge (see Section 5). Secondly, that the outflow rate M is likely to vary 
with variable geometry of irradiation, i.e. with ^d,i^ or with (|)*^^. 

Therefore the shape of the light curve can be formally written as 

l{^nsh) ~ M(^:,J \fe{<^nsh) + ] , (29) 

where (|)*^,/^ is the phase at the moment of irradiation, i^nsh - the phase at the moment 
of impact, and "x" in front of ft is intended to represent the higher value of the 
impact parameter Ay-^ in the case of the overflowing parts of the stream. 

Limiting our discussion to qualitative considerations we recall (see Fig. 3 in 
Section 5) that variations of show one maximum (and one minimum) per cycle 
while those of - two maxima (and two minima) per cycle. Variations of M are 
also expected to show two maxima (and two minima) per cycle. Taking this into 
account we conclude that (1) the contribution from overflowing parts of the stream 
ft must be dominant (due obviously to x > 1 ; see above), and (2) the contributions 
from variations of fe and M make the shape of the light curve significantly different 
from a simple cosine wave. 

The second prediction is confirmed (at least qualitatively) by the observed 
shapes of negative superhump light curves which in nearly all cases (Barett et al. 
1988, Fig.4; Patterson 2001, Fig.3; Patterson et al. 1997, Fig.4; Patterson et al. 
2002, Figs.8,10,15; Stanishev et al. 2002, Fig.3) are clearly non-cosinusoidal. 

73. The Light Variations with Pprec 

The negative superhump maximum, which occurs - by definition - at <^„sh = 0, 
is produced when the stream is overflowing the top part of the disk at the position 
angle 6 = 0. Using Eqs.(19) and (20) we get 

do = a. (30) 

As mentioned in the Introduction, the mean luminosity of the disk varies with 
Pprec- The maximum of those variations occurs when Qd,obs = 0- Combining 
Eqs.(15) and (30) we then predict that it should ocur at 
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^prec = a ~ 0.014 . (31) 

This agrees nicely with observations (e.g. TT Ari - Semeniuk et al. 1987, TV Col 
- HelUer 1993 and PX And - Stanishev et al. 2002) which show that the maximum 

occurs at ^prec ~ 0. 

8. Discussion 

The model presented in this paper appears quite simple and self-consistent. 
Its predictions compare favorably with observations. There are several problems 
and questions, however, which should be answered prior to considering it as fully 
acceptable. 

(1) What makes the originally coplanar disk to become tilted? The qualitative 
answer to this question may be quite simple: when some part of the outer disk - 
due to a random fluctuation - is deflected above or below the orbital plane, then the 
mechanism described above can begin to operate. What conditions, however, are 
necessary for producing the positive feed-back (especially for condition given by 
Eq.25 to be fulfilled)? 

(2) What is the cause of transitions from a tilted disk (with negative super- 
humps) to a coplanar disk (with common superhumps) - and vice-versa - observed 
in many systems (e.g. in V603 Aql or TT Ari; see references to Table 1)? 

(3) How can we explain the simultaneous presence of negative and common 
superhumps in some systems (e.g. in V603 Aql or V503 Cyg; see references to 
Table 1)? And why are such cases so rare? 
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